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VARIABLE RATE CONTINUOUS MODE SATELLITE MODEM 



FIELD OF THE INVENTION 



The present invention relates generally to communications and more particularly 
relates to a variable rate continuous mode modem suitable for satellite communications. 



BACKGROUND OF THE INVENTION 



Recently there has been a widespread proliferation of satellite communications 
systems throughout the world. A large number of consumers have acquired satellite 
communications equipment for home reception of satellite broadcast television signals in 
addition to many businesses having acquired such equipment for the purpose of receiving 
and/or transmitting useful business information. 

Typically, satellite communications systems have been constructed using a 
combination of RF circuitry, IF and baseband circuitry. A common trend in the art is to 
implement the IF and baseband modules using one or more FPGAs and/or ASICs. A 
disadvantage of this is that it is difficult or even impossible to perform some techniques that 
are only possible to perform in software. Attempting to perform these techniques in circuitry, 
if possible, would result in large, complex and expensive silicon. 

Thus it is desirable to have a satellite modem wherein the receiver portion is mainly 
implemented in software so as to reduce cost, making necessary changes much simpler to 
perform and potentially reducing the size of the receiver. 
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SUMMARY OF THE INVENTION 
The present invention overcomes the disadvantages of the prior art by providing a 
satellite modem whose baseband modules can be implemented mostly in software running on 
a suitable processor such as a high speed digital signal processor (DSP) which are commonly 
5 available today. The modem includes a RF portion and a baseband portion. The digitized 
output of a matched filter is input to the baseband processing portion that consists of software 
executing on a DSP. 

The modem comprises an antenna coupled to an upconverter/downconverter. The 
upconverter/downconverter is coupled to a RF transceiver which functions to receive a 

10 transmit signal from a baseband transmitter module and generate a signal for input to the 
upconverter/downconverter. The RF transceiver also functions to receive a downconverted 
signal and output an IF signal for input to an IF module. The IF functions to generate I and Q 
outputs that are input to a digital baseband module. The baseband module can be 
implemented as software executing on a DSP. 

15 A novel aspect of the present invention is the method of performing phase acquisition. 

In particular, the invention disclosed a technique of modulation wipe off wherein the effects 
of modulation are removed. The technique involves rotating all the symbol vectors into a 
single quadrant. This is done in order to avoid the case where the constellation points are 
spread over more than one quadrant due to excessive noise being present in the channel. 

20 Another novel aspect of the invention is the signal detection technique. The signal 

detection technique performs signal detection and frequency acquisition in the presence of 
multiple signals in accordance with a method described hereinbelow. 

Still another novel aspect of the present invention is the timing acquisition technique 
that includes an unwrapping technique. This functions to correct timing errors before 

25 averaging is performed. In addition, linear changes to the time estimates are performed using 
LSR techniques. 

The satellite modem of the present invention, wherein the baseband receiver and 
transmitter are mainly implemented in software, provides for a reduction in cost, makes 
performing necessary changes much simpler and reduces the size of the modem. In addition, 
30 the invention represents a significant step toward the goal of a complete software radio. 

There is thus provided in accordance with the present invention a satellite modem, 
comprising a housing, a antenna for receiving and transmitting radio frequency (RF) signals, 
an upconverter/downconverter coupled to the antenna and adapted to upconvert RF signals 
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from a first frequency band to a second frequency band, and adapted to downconvert RF 
signals from the second frequency band to the first frequency band, an RF transceiver 
coupled to the upconverter/downconverter and adapted to receive an output transmit signal 
and to generate an output RF signal therefrom, the RF transceiver adapted to receive an input 
5 RF signal from the upconverter/downconverter and generate an input receive signal 
therefrom, an intermediate frequency (IF) module adapted to receive the input receive signal 
and generate I and Q signals in response thereto, a baseband module adapted to receive the I 
and Q signals and to generate receive data in accordance therewith and a baseband module 
adapted to generate the output transmit signal in accordance with a transmit data signal input 
10 thereto. 

There is also provided in accordance with the present invention a receiver baseband 
apparatus, comprising input means adapted to receive an I and Q signal, an I matched filter 
adapted to receive the I signal and generate an I filtered output therefrom, a Q matched filter 
adapted to receive the Q signal and generate a Q filtered output therefrom, a processor 

15 programmed to perform automatic gain control (AGC) and generate an AG control signal 
therefrom, perform timing detection and generate an A/D clock control signal therefrom, 
perform phase detection and generate a voltage controlled oscillator (VCO) control signal 
therefrom, a decoder adapted to receive the I output signal and the Q output signal from the 
processor and to generate a decoded output therefrom, a deinterleaver adapted to generate a 

20 deinterleaved output in accordance with the decoded output signal input thereto, a forward 
error correction decoder adapted to generate output receive data in accordance with the 
deinterleaved output signal input thereto and a controller adapted to manage and control the 
input means, I matched filter, Q matched filter, the processor, the decoder, the deinterleaver 
and the forward error correction decoder. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
The invention is herein described, by way of example only, with reference to the 
accompanying drawings, wherein: 

Fig. 1 is a block diagram illustrating the RF, IF and baseband portions of the variable 
5 rate satellite modem of the present invention; 

Fig. 2 is a block diagram illustrating the baseband module portion of the variable rate 
satellite modem in more detail; 

Fig. 3 is a flow diagram illustrating the main software processing loop of the variable 
rate modem; 

10 Fig. 4 is a flow diagram illustrating the processing performed upon power up or after 

the communication link is broken; 

Fig. 5 is a flow diagram illustrating the initial automatic gain control processing 
portion of the variable rate modem; 

Fig. 6 is a chart illustrating the results of correlation between the spectral estimation 
15 , and the spectral mask of the expected signal with both desired and adjacent signals having bit 
rates equal to 9600 bps; 

Fig. 7 is a chart illustrating the correlation result with the desired signal having a bit 
rate equal to 9600 bps and an interfering signal having a bit rate of 128 Kbps; 

Fig. 8 is a flow diagram illustrating the signal detect and frequency acquisition 
20 processing portion performed upon power up; 

Fig. 9 is a flow diagram illustrating the timing, AGC, frequency and phase acquisition 
processing portion of the variable rate modem; 

Fig. 10 is a flow diagram illustrating the fine frequency estimation processing portion 
of the variable rate modem; 
25 Fig. 1 1 is a flow diagram illustrating the coarse phase acquisition processing portion 

of the variable rate modem; 

Fig. 12 is a flow diagram illustrating the fine phase acquisition processing portion of 
the variable rate modem; 

Fig. 13 is a block diagram illustrating tracking portion of the variable rate modem; . 
30 Fig. 14 is a flow diagram illustrating the pre-track processing portion of the variable 

rate modem; 

Fig. 15 is a flow diagram illustrating the AGC tracking portion of the variable rate 
modem; 
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Fig, 16 is a block diagram illustrating a model of the timing tracking portion of the 
variable rate modem; 

Fig. 17 is a block diagram illustrating the phase tracking portion of the variable rate 
modem; and 

Fig. 18 is a flow diagram illustrating the tracking processing portion of the variable 
rate modem. 
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DETAILED DESCRIPTION OF THE INVENTION 
Notation Used Throughout 



The following notation is used throughout this document. 



lerm 


Definition 


AGC 


Automatic Gain Control 


ASIC 


Application Specific Integrate Circuit 


BER 


Bit Error Rate 


DDS 


Direct Digital Synthesizer 


DEC 


Circular Dnft Function of EC 


DFT 


Discrete Frequency Transform 


DSP 


Digital Signal, Processor 


ES 


Estimated Spectrum 


. T7T7 


Frequency Estimate 


FFT 


Fast Fourner Transform 


FIR 


Fimte Impulse Response 


FPGA 


Field Programmable Gate Array 


lr 


Intermediate Frequency 


IIR 


Infinite Impulse Response 


LIM 


Limit of Uncertainty Range 


LPF 


Low Pass Filter 


LSR 


Least Squares Regression 


NCO 


Numerically Controlled Oscillator 


NFFT 


Number of FFT points 


PI 


Peak Indexes 


PLL 


Phase Lock Loop 


QPSK 


Quadrature Phase Shift Keying 


RF 


Radio Frequency 


RMS 


Root Mean Square 


SNR 


Signal to Noise Ratio 


TS 


Theoretic Spectrum Mask 


VCO 


Voltage Controlled Oscillator 



5 General Description 

The present invention is a variable rate continuous mode satellite modem capable of 
transmitting and receiving at a symbol rate in the range of 9.6 K to 256 K symbols/sec. Note 
that the bit rate depends on the coding rate. For example, at a code rate of V2, with QPSK, the 
bit rate is equal to the symbol rate. At a code rate of with QPSK, the bit rate is 50% more 
10 than the symbol rate. To aid in understanding the principles of the present invention, the 
description of the processing performed in the modem is divided into five main portions: (1) 
the main loop including a state machine for defining acquisition, tracking, etc., (2) signal 
detection and initial frequency estimation, (3) acquisition of timing, frequency and phase, (4) 
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pre-tracking (AGC, timing and phase tracking) with a phase loop having a larger bandwidth 
and phase ambiguity to be solved and (5) tracking of the AGC, timing and phase. Additional 
processing includes the soft decision data detection. 

The main loop will now be described in more detail. A block diagram illustrating the 
5 RF, IF and baseband portions of the variable rate satellite modem of the present invention is 
shown in Figure 1. The modem, generally referenced 10, comprises a housing (not shown), 
an upconverter/downconverter 13 coupled to an external antenna 12 and a RF transceiver 14 
coupled to the upconverter/downconverter 13. The upconverter functions to convert the 
output of the RF transceiver 14, e.g., L band, to the RF output frequency, e.g., C or Ku band. 

10 The downconverter performs downconversion of the signal received over the antenna from 
RF frequencies, e.g., C or Ku band, to L band for input to the RF transceiver 14. The RF 
transceiver downconverts the received RF signal to an IF signal. The IF signal output of the 
RF transceiver 14 is input to the IF module 16 which functions to generate the I and Q 
baseband signals. The baseband module 18 demodulates the I and Q signals to yield the 

15 receive Rx data! The baseband module 18 performs signal processing tasks on the I and Q 
data from the IF module 16 such as Viterbi decoding, de-interleaving and Reed Solomon 
decoding 

In the transmit direction, the baseband transmitter 20 generates a baseband transmit 
signal in response to the input transmit Tx data. The baseband transmit signal is input to the 
20 RF transceiver 14 which upconverts the signal to an RF frequency signal which is then input 
to the external antenna 12. 

Baseband Module 

A block diagram illustrating the baseband module portion of the variable rate satellite 
modem in more detail is shown in Figure 2. The baseband module 18 receives the baseband I 

25 and Q signals. The analog I signal is converted to digital by A/D converter 21. The 8-bit 
digitized I signal is input to a matched filter 23. The analog Q signal is converted to digital 
by A/D converter 22. The 8 -bit digitized Q signal is input to a matched filter 24. The 
matched filters 23, 24 are implemented in hardware but may alternatively be implemented in 
software. The I and Q outputs of the matched filters 23, 24 are input to processor 25. 

30 In one embodiment of the invention, for example, the RF receive frequency for C 

band operation is approximately 4 GHz and the RF transmit frequency is approximately 6 
GHz. The RF transceiver 14 do wnconverts/upcon verts to/from an IF frequency of 
approximately 900 MHz. The IF module 16 downconverts the IF signal to analog I and Q 
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baseband signals within the bandwidth range of approximately 0 to 170 KHz. The baseband 
module 1 8 demodulates the baseband I and Q signal to a soft decision digital data signal. 
Note that the IF frequency, between the RF transceiver 14 and the IF module 16, may be 70 
MHz ±18 MHz or 140 MHz ±36 MHz. The RF transceiver 14 may comprise a 
5 commercially available transceiver suitable for operation with the satellite modem of the 
present invention. For Ku band operation, the RF receive frequency is approximately 12 
GHz and the RF transmit frequency is approximately 14 GHz. The RF transceiver 
downconverts/upconverts to/from an IF frequency within the range of approximately 950 to 
1450 MHz. 

10 In an alternative embodiment, the signal from the RF transceiver constitutes a first IF 

which is input to a modified IF module that outputs a second IF signal. This second IF signal 
is converted to digital via an A/D converter and input to a digital converter that functions to 
generate I and Q signals. 

The processor 25 performs the majority of the signal processing functions in the 

15 modem 10 such as activity detection, timing acquisition, frequency acquisition, phase 
acquisition and tracking. The processor 25 functions to output 3-bit I and Q signals that are 
input to a Viterbi decoder 26. The output of the Viterbi decoder 26 is input to a deinterleaver 
27. The output of the deinterleaver 27 is input to a Reed Solomon decoder 28 which 
functions to output the receive Rx data. 

20 A controller 29 controls and manages the operation of the baseband module and clock 

circuitry 3 1 functions to provide the various clocks needed by the elements of the baseband 
module 18. 

Main Processing Loop 

A flow diagram illustrating the main software processing loop of the variable rate 
25 modem is shown in Figure 3. Upon power up (step 32) a series of start up tasks are 
performed. These include activity detection (step 34), acquisition (step 36) and pre-tracking 
(step 37). The activity detection (step 34) performs signal detection, AGC and coarse 
frequency acquisition. In this step, the modem is listening for a signal from the transmitter. 
Hie frequency /tuning offset is reduced to a suitable level such that when the PLL is started, 
30 the errors are relatively small and the PLL is capable of locking onto the signal. During the 
activity detect phase, no matched filtering is performed. In addition, the frequency is variable 
with an approximately frequency offset of 30 KHz at a low rate and with a bandwidth of 
approximately 7 to 8 KHz. 
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During acquisition (step 36), AGC is performed as well as timing, fine frequency and 
phase acquisition. Next, during the pre-tracking phase (step 37), phase tracking acquisition is 
performed wherein a phase estimate is performed every four symbols. 

Once the start up phase is completed, tracking is performed (step 38). This step 
5 includes AGC, timing and phase tracking. The above steps are repeated when tracking is lost 
on the signal or after each time the modem is powered up. 

Activity Detection 

Each of the above steps is described in more detail hereinbelow beginning with the 
activity detection phase. A flow diagram illustrating the activity detect method (step 34 
10 Figure 3) is shown in Figure 4. The following steps are performed upon power up or after a 
catastrophic event whereby the communication link is completely broken. The initial 
automatic gain control (AGC) step includes measuring the energy at the output of the A/D 
converters 21,22 (Figure 2), i.e., at the input buffer (step 40). 

Next, the decimation of the input signal is performed (step 42). Decimation by a 
15 factor of 1, 2, 4 or 8 is performed. Signal detection and initial coarse frequency acquisition 
and estimation is performed (step 44). If signal activity is not detected, the method loops 
back and repeats steps 40, 42, 44 until a signal is detected (step 46). 

Each of the steps of the activity detection method will now be described in more 
detail beginning with the initial AGC step 40. A flow diagram illustrating the initial 
20 automatic gain control (AGC) processing portion of the variable rate modem is shown in 
Figure 5. First, N samples of the input signal are collected whereby N may be 64, 128 or 256 
(step 50). The unfiltered samples are expressed as 

z n =I n +JQ n (1) 
The mean square energy Ems is then calculated using the following equation (step 5 1). 

25- ^^EN ( 2 ) . 

where G Q is the gain of the AGC and IG G is the inverse gain. The inverse gain IG C is set equal 
to 256 * E ms for scaling the RMS of the A/D converter input at 1/16 maximum input voltage 
level (step 52). The value of G 0 is determined from a lookup table (step 54). The G Q value is 
then written to the AGC hardware (step 56). Signal detection is then performed (step 58). 
30 Signal decimation (step 42 Figure 4) is performed before the actual processing in . 

order to improve the frequency resolution at lower data rates. The resolution achieved is a 
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particular fraction of the symbol rate. It is preferable to have between 4 and 10 samples per 
symbol, except for the case of very low data rates. This means that the signal is relatively 
wideband in the normalized frequency domain thereby enabling a reasonable resolution 
without requiring very long FFTs. The resolution can be expressed as 

binwidth _ f s I N 



Symbol Rate R b 



nss 



N 



(3) 



where 

f s is the sampling rate 

Rb is the symbol rate 

nss is the number of samples per symbol 

N is the number of points in the FFT 
The parameters of the decimation for activity detection, also referred to as gear shifting, is 
presented below in Table 1. 



Table 1 



Symbol Rate 


Sample Rate 


Max BW 


Decimation 


Number of 


Normalized 


Relative 


(KHz) 


at A/D (KHz) 


(KHz) 


Factor 


samples per . 


Cutoff 


Resolution (%) 










symbol 


Frequency 




9.6-12 


768 


38.1 


8 


8 - 10 


0.0496 


3.1 -3.9 


12-16 


864 


41 


8 


6.75 - 9 


0.0475 


2.63-3.5 


16-24 


928 


46.2 


8 


4.83 - 7.25 


0.04978 


1.88-2.8 


24-32 


512 


51.6 


4 


4-5.33 


0.1008 


1.5-2 


32-64 


512-1024 


51.6 


4 


4 


0.1008 


1.5 


64 - 128 


512-1024 


73.2 


2 


4 


0.1430 


1.5 


128-256 


512-1024 


116.4 


1 


4 


0.2273 


1.5 



15 

The maximum bandwidth is determined by the bandwidth of the higher rate plus a 
maximum frequency offset of 30 KHz. 

The relative required resolution is determined from the limitations of the phase 
acquisition algorithm. The phase acquisition algorithm operates on blocks of eight symbols, 
20 trying four hypotheses for the axis, each rotated by 22.5 degrees from one another. 

Therefore, over eight symbols, the maximum frequency offset causes less than 10 to 
20% of 22.5 degrees, i.e., less than 2.5 to 5 degrees phase shift. This implies a resolution of 
much less than one degree per symbol. This corresponds to a resolution of less than 0.3% of 
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the symbol period. This resolution, however, is very difficult to obtain, thus a fine frequency 
acquisition is performed prior to the phase acquisition. This yields resolutions between 1.5 to 
3% of the symbol period. 

The decimation is performed by filtering utilizing half band filters, which results in 
5 decimation by a factor of two. In addition, these filters are relatively inexpensive to 
implement since (1) approximately half the coefficients are zero, (2) the filter order is 
relatively low due to a wide transition band between the pass and stop band and (3) only 
every second output need be computed due to the decimation by two. The design of half 
band filters is described in more detail in "Multirate Systems and Filter Banks," P.P. 
10 Vaidyanathan, Prentice Hall 1993; "A Trick for the Design of FIR Half Band Filters," P.P. 
Vaidyanathan and T. Nguyen, IEEE Transactions on Circuits and Systems, Vol. CAS-34, pp. 
378-389, April 1 989, both of which are incorporated herein by reference. 

In the embodiment presented herein, a filter having a length of 9 for the first stage is 
sufficient, while a filter length between 11 and 13 is sufficient for the second stage (or the 
15 single stage when decimation by 2 is performed). 

The coefficients for the filters are presented below in Tables 2 through 4. 
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Table 2: Coefficients For Decimation By Factor Of 8 Using Three Filters 



-0.0336 


0 


Filter #1 

0:2835 0.5000 0.2835 


0 


-0.0336 


-0.0422 


0 


Filter #2 

0.2902 0.5000 0.2902 


0 


-0.0422 


0:0131 


0 


Filter #3: Columns 1 through 7 
-0.0236 0 0.0463 


0 


-0.0949 


0 


0.3144 


Filter #3: Columns 8 through 14 
0:5000 0.3144 0 


-0.0949 


0 


0.0463 


o 


Filter #3: Columns 15 through 19 
-0.0236 0 0.0131 






Table 3: Coefficients For Decimation By Factor Of 4 Using Two Filters 


-0.0426 


0 


Filter #1 

0.2904 0.5000 0.2904 


0 


-0.0426 


0.0146 


0 


Filter #2: Columns 1 through 7 
-0.0246 0 0.0472 


0 


-0.0956 


0 


0.3147 


Filter #2: Columns 8 through 14 
0.5000 0,3147 0 


-0.0956 


0 


0.0472 


0 


Filter #2: Columns 15 through 19 
-0.0246 0 0.0146 






Table 4: Coefficients For Decimation By Factor Of 2 Using One Filter 


0.0169 


0 


Filter #3 : Columns 1 through 7 
-0.0691 0 0.3041 


0.5000 


0.3041 


0 


-0.0691 


Filter #3 : Columns 8 through 1 1 
0 0.0169 
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Signal Detection and Frequency Acquisition 

This section describes the method used by the variable rate modem for signal 
detection and frequency acquisition in the presence of multiple signals in the spectrum. 
According to Intelsat standards, two signals can be placed at frequencies that are a distance of 
5 G.7 * (symbol_ratei + symbol_rate2) as measured from center frequency to center frequency 
E 

and whereby the — — levels differ by at most 7 dB. 

The methods of the present invention utilize FFT periodograms for spectral estimation 
followed by cyclic correlation with the spectral mask of the expected signal in order to detect 
the center of the signal. Figures 6 and 7 are charts illustrating the results of correlation 
10 between the spectral estimation and the spectral mask of the expected signal. The result of the 
correlation for the Intelsat scenario, for example, is shown in Figures 6 and 7. In both 
figures, the expected signal has a power level lower by 7 dB than adjacent signals and has a 
data rate of 9600 bps. 

With reference to Figures 6 and 7, multiple peaks can be detected in the correlation. 
1 5 In contrast to the single signal scenario, the maximal peak is not the sole candidate for the 
central frequency of the expected signal. Preferably, peaks in the correlation should be 
examined to determine the potential locations of the desired signal. 

Spectrum Estimation 

A flow diagram illustrating the signal detect and frequency acquisition processing 
20 portion (step 44 Figure 4) performed upon power up is shown in Figure 8. The first step 
performed is to calculate the spectrum of the received signal (step 60). The estimated 
spectrum (ES) is represented as an array or size 256. The following pseudo code describes 
this step. 

Spectral Estimation Pseudo Code 



25 



1. reset the ES array to zero; 

E 

2. for j=0 to num_reps - 1 (num_reps = 32 for low — — levels and 16 for regular operation) 



get z , a vector of 256 received signal samples, z i — I i + jQ n i = 0..255 
ES = ES + \FFT(ztf 



30 endfor 

ES 



3. calculate ES = 



num _ reps 
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4. FFT shift ES: 

(ES(i- 128), / = 128*255 .„ , 

ES'U) = \ , where ES' is the FFT shifted version of 

W [£:5(i + 128),/ = 0*127 

£S, throughout the remainder of this section, however, ES' is denoted as ES; 

5. remove DC using: 

£S(128) = SiM ; 



Signal Presence Determination 
Next, the presence of a signal is determined using the following pseudo code (step 

62). 

10 Signal Presence Determination Pseudo Code 

1. smooth the Estimated Spectrum using an 11 coefficient low pass filter (LPF) 
presented below in Table 5; 

15 ESS(0=^ES(i- j)hU),i = 0*255 

where ES(i) = 0 for /< 0 and i> 255; 

Table 5 : Coefficients For Low Pass Filter 
h_ 5 bu h_3 h_2 h-i ho hi h 2 h 3 lu h 5 

0.010 0.0249 0.0668 0:i249 0.1756 0.1957 0.1756 0.1249 0.0668 0.0249 0.010 

20 2. calculate RES, the resolution of one bin in the FFT 

_ SamplingFrequency 

KEo = — 

NFFT 

where NFFT is defined as the number of FFT points; 

3. calculate chk_rng, the range of indexes to be checked in the FFT 

f 4 X 

i max(3xl0 ,0.7 -symbol rate 
chk __rng ; = ceil\ • 

V y 

25 . 4. calculate max_chk and min_chk 

max_chk = max{£S 5 (z)C^"f ng 

min_chk = minfe (0}!!^wT^ 
5. calculate act ratio 
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max chk 

act _ ratio = — ; — 

mm_chk 

If the value of act_ratio is greater than 1.75, it is determined that a signal is present (step 64). 
Otherwise it is determined that no signal is present. 

5 Correlation Calculation 

The next step is to perform a correlation calculation (step 66). The following 
definitions apply: 

ES: estimated spectrum 
TS: theoretic spectrum mask 
10 EC: an array of size 256 representing the correlation function between the non 

smoothed ES and the TS 

LIM: the limit of the range of uncertainty 

The correlation is determined using the following pseudo code (step 66). 
15 Correlation Determination Pseudo Code 

1. calculate LIM using the following equation; 

/ 3 0,000 + 0.125 symbol rate^\ 

LIM = ceiA — - = 

I, RES J 

2. calculate the following; 

20 EC^i- len ^KTS) Jj = /CTg g^( ( ^ + i) mod256 yrS(k)MS - Lim < i < 128 + Lim 

Peak Determination 

The next step is to perform a peak determination algorithm (step 68). The following 

definitions apply: 

25 DEC: circular difference function of EC 

PI: list of peak indexes 

N WIN: the minimal distance between peaks 



30 



The peak determination is determined using the following pseudo code. 



Peak Determination Pseudo Code 



1 . calculate DEC using the following equation; 

z, £C( o4 £C(0) - £C(255),i=0 

W \EC{i)-EC{i-\),i = \..255 
35 2. find PI, for the list of indexes i which fulfill the conditions / e PI 
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(DEC(i) > 0 and DEC(i + 1) < 0, i = 0. .254 
\dEC(255) > 0 and DEC{6) < 0, / = 255 

3. find the value of the minimum peak and eliminate from the list those peaks that are 
10 dB (0. 1) smaller than the maximum peak; 

4. calculate N_WIN using the following; 

NFFT^ 



N WIN = < 



floor] 0.35 



nss 



symbol _ rate < 64 KHz 



floor 



NFFT^ 



0. 1 75 symbol _ rate > 64 KHz 



V nss J 

5. If two or more peaks are closer than NJWTN, keep the highest between them in the 
list. The order of the removal is from the highest peak to the lowest. Note that it is well 
known that the center frequencies of two signals are separated by the sum of half the signal's 
bandwidth, i.e., 0.7 {symbol _rate } ^symbol _rate 2 ). In practice, since (1) the peaks of 
10 interferers are not accurately detected, (2) the peak of a small rate signal tends to get closer to 
the center of the high rate signal, and (3) the rate of the interferer is not known, this range 
should be decreased to the value given above. 

Symbol Rate Calculation 
The next step is to calculate the symbol rate (step 70). The following definitions 

15 apply: 

BW2: one half of the bandwidth of the signal 

The symbol rate is calculated utilizing the following pseudo code. 

20 Symbol Rate Calculation Pseudo Code 

1. calculate BW2 using the following; 

f 0.675 - symbol rate\ 

BW2 , cal > rks — j 

2. calculate the following; 

25 No _ level ^ = Afea«({^5'(/)|0.85 -min_cM < ES(i) < 1.25 -min_cM,0 < i < 255) 

3. for /= 1 to the length of PI: 

3.1 define Area as: 

Area = {/|P/(/) — BW2 <i < PI (I) + BW2 

3.2 define Area2 as: 



30 Area = { i 



PI(!)--BW2<i<PI(l) + -BW2 
8 8 



3.3 find max val: 



max val = max ES s (k) 

k&Areal 
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3.4 redefine the indexes of Area2 as k such that ES s (k) < maX Va * 



2 

f ^ 1 

Area2 



PI{1) - 1 BW2 < i < PI(l) + \bW2, ES s (i) < maX Vd 



8 8 * 2 

3.5 calculate mean_val using; 

mean_val = 

number of indexes in Area2 

5 3.6 find Fu, the smallest index in Area greater than PI(1), such that 

ES S (F u ) < 0.5 • (mean _ val + No _ /eve/) 

3.7 find Fl, the highest index in Area smaller than PI(l), such that 
ES s (F t ) < 0.5 • (mean_val + No _ level) 
10 3.8 if both Fu and Fl are found then calculate the estimated symbol rate 

(es_symbol_rate) corresponding to PI(1) using the following: 
es _ symbol _ rate ~ (F u — Ff} RES 

4. end for 

5. remove from the list the peaks having es_symbol_rate higher than 1.2 times the 
15 expected symbol rate or lower than 0.75 times the expected symbol rate. 

If the list of peaks is exhausted, the method continues with step 60 (step 72). If the 
list of peaks is not exhausted, the peak closest to the center is selected (step 74). 

Frequency Acquisition 

20 The next step is to perform frequency acquisition (step 76). The following definitions 

apply: 

FE: the normal frequency estimation 

pi; the peak index where a candidate signal is centered 

25 The frequency acquisition is performed utilizing the following pseudo code. 

Frequency Acquisition Pseudo Code 

1. the normal frequency estimation (FE) is calculated using the following 
30 FE = (pi — 128) - RES 

2. FE is output to the synthesizer and the central frequency is shifted in order that the 
frequency shift be compensated. 

Communications is then attempted (step 78). If communications are not successful, 
35 the particular peak is deleted from the list and the method continues with step 72). If 
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communication is successful, the method terminates. Timing acquisition and fine frequency 
acquisition are then performed (described in more detail below). 

Timing, Frequency and Phase Acquisition 

A flow diagram illustrating the timing, AGC, frequency and phase acquisition 
5 processing portion (step 36 Figure 3) of the variable rate modem is shown in Figure 9. A 
matched filter begins operation at a normal gear shift at 4 samples per symbol (step 90). The 
matched filter is implemented in hardware and operate without a phase rotation. 
Alternatively, the matched filter can be implemented in software. The matched filter 
hardware controls the clocking of the A/D converters. 
10 The normal gear shift ratios are presented below in Table 6. 



Table 6: Normal Gear Shift Table 



Index 


Low Symbol 
Rate (KHz) 


High 

Symbol Rate 
(KHz) 


Samples/ 
Symbol 


Low 

Frequency 
(KHz) 


High 

Frequency 
(KHz) 


7 , 


9.6 


16 


64 


614.4 


1024 


6 


16 


32 


32 


512 


1024 


5 


32 


42 


24 


768 


1008 


4 


43 


64 


16 


672 


1024 


3 


64 


85 


12 


768 


1020 


2 


85 


238 


8 


680 


1024 


1 


128 


170 


6 


768 


1020 


0 


170 


256 


4 


680 


1024 



First, 128 samples are collected. At 4 samples per symbol, this is equivalent to 32 
15 symbols. The first AGC acquisition is performed on the block of 128 samples (step 92). 
During this step, the following is calculated. 

1 n 2 

where N = 128 samples at 4 samples per symbol and the AGC value is calculated 
utilizing the following equation: 

Ems 



20 AGC = 



2xl6 2 



The next step is to perform timing acquisition (step 94). The following definitions 

apply: 

/block of symbols': 16 contiguous symbols 
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ti: DFT estimate of a block of symbols, i - 1 ..N 
N: number of DFT estimates 

'group of DFT estimates': group of N DFT estimates 
Tj: average estimate of group of N DFT estimates 
5 K: number of average estimates Ti, which is the number of points used as input to the 

Least Squares Fit 

The timing acquisition is performed utilizing the following pseudo code. 
In connection with the following Pseudo Code, the following assumptions are made. 
10 (1) blocks of 16 contiguous symbols at 4 samples per symbol are available, i.e., 16 * 4 = 64 
continuous samples; (2) the next block of 16 symbols arrive after a break of 16 * 3 = 48 
symbols; and (3) the maximum drift between the receiver and the transmitter clocks is 
approximately 100 ppm. 

The method of the invention utilizes the timing recovery technique described in M. 
15 Oerder and H.Meyr, "Digital Filter and Square Timing Recovery," IEEE Transactions on 
Communications, Vol. COM-36, pages 605 to 612, May 1988, incorporated herein by 
reference. The method of the invention is given briefly below. 

1 . Perform steps 2 through 5 for K groups made up of N DFT estimates, each 
calculated from blocks of 16 symbols. 
20 2. Calculate timing estimate t u based on the DFT technique for 16 contiguous 

symbols, i = 1,..,N, thereby obtaining N DFT estimates each based on 16 symbols. 

3. Make a histogram (3 bins in the range of -0.5, 0.5) of the N DFT estimates tj 
obtained in step 2. 

4. Classify the timing range (0, %, l A, -V4) the group the N estimates are in. The 
25 decision is derived from logic based on the histogram in step 3. 

5. Unwrap N DFT estimates based on a range decision derived from logic and 
calculate their average T;. 

6. Unwrap the K average estimates T; and do a least square fit to the K averages to 
produce the final estimate. 

30 

Timing Acquisition Pseudo Code (in more detain 

1 . for k = 1 to K perform steps 2 through 5 inclusive )groups of average estimates Ti); 

2. calculate the timing estimate ti based on a DFT algorithm for 16 contiguous 
35 symbols, i = 1..N, thereby obtaining N DFT estimates each based on 16 symbols contiguous 

at 4 samples per symbol with no discontinuity in time. 

for i= 1 toN 
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e r = jr{log[l + SA^|z^ 

^- = ^ ^ log[l H- • |2:(4 y + H- 16 * 4 • (/ - l))f ]-h log[l + .Wi? • |z(4 y + 3 + 16 - 4 - (/ - 1))| 2 ] 



endfor (i = 1 to N) 
where 

E 



SNR = 



s_ 
"a 



and E s represents energy per symbol. 

The preliminary estimate is then given by the following 

i — l..N DFT estimates 



1 . *i 
t. = arctan — 

2; pi Uv 



Note that the arctan function can be approximated using the least squares method and 
the logarithm is the natural logarithm to the base 'e'. 
10 3. construct a histogram of the N DFT estimates ti wherein 3 bins in the range -0.5 to 

0.5. 

3.1 count the number of estimates t; in each bin as defined by 
bin 1 = [-0.5, -1/6] 

bin 2 = [-1/6, 1/6] 
15 bin 3 = [1/6, 0.5] 

3.2 sort counts of the three bins in descending order and save their respective indices 
as follows: 

max = value of the maximum count 
20 mid = value of the middle count 

min = value of the minimum count 

Imax = index of the bin with the maximum count 
Imid - index of the bin with the middle count 
25 I m in — index of the bin with the minimum count 

4. classify the portion of the timing range the group of N estimates are in; 

range 1 = around zero, i.e., [-1/8, 1/8] 
30 range 2 = around !/4, i.e., [1/8, 3/8] 

range 3 = around X A or -14, i.e., [- l / 2 , 3/8] and [3/8, X A\ 
range 4 = around - l A, i.e., [3/8, -1/8] 

if (max > 10) & (mid + min < 6) then 
35 if I max = 2 then 

range =1; time round 0 
elseif I max = 1 

range = 4; time negative 

else 
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else 



10 end 



range = 2; time positive 

end 

if Imin = 2 & (max - mid < 4) 
range = 3; time round Vi 

elseif I min = 1 

range = 2; time positive 

else 

range = 4; time negative 



5, unwrap N DFT estimates based on the range decision and calculate their average T k 
using the following pseudo code; 

15 Average TV Calculation Pseudo Code 

5.1. initialize thresholds: 
thrmid = 0.25 
thrhlf=0.05 
20 5.2. calculate 

mean = > % average of DFT estimates before unwrapping 

5.3. do for group of estimates i = 1 to N 

sgn = 0 % case on value of the variable 'range' 
case 1 

25 NOP 
case 2 

if ti < rthrmid then sgn = 1 ; 
case 3 
if mean > 0 

30 if ti <-thrhlf then sgn = 1; 

else 

if ti > thrhlf then sgn = -1; 
end 
case 4 

35 if ti > thrmid then sgn = -1 ; 

end; 

end do for group of estimates 

5.4. calculate the average estimate of the group of N DFT estimates as follows 



40 T. - , f = 1 . JV , % calculate average estimate 

N 

endfor (k = 1 to K) 

6. unwrap the K average estimates T t and do a least square fit to the K averages to 
yield the final estimate; 

45 We assume that Ti depends linearly on time. The observation interval, however, is 

limited to the range [-0.5, 0.5]. 
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6. 1 unwrapping the K average estimates 

The first step is to unwrap the K average estimates and receive new values Tj for 

K-l K + \ ~ ■ „ 

i= foroddK. 

2 2 

5 

Unwrapping Pseudo Code 

1. Detect positive and negative jumps in timing that are larger than 0.5. A positive 
jump is detected if 
10 TO> + l)-r(;»>0.5 

A negative jump is detected if 

ro« + i)-rc/n)<o.5 

where T(j) are the time estimates and j is the time index. 

15 2. For positive jumps at index jp, subtract 1 from all subsequent time estimates (with 

indices > jp. 



20 



30 



3. For negative jumps at index jn, add 1 to all subsequent time estimates (with indices 

> jn. 

6.2 Calculate estimates of the slope a and the average b of the Tj, i = 1 to K, using a 
least square fit estimation 



The estimate is calculated to be: 
25 ' T = a < K - l Kb 



with no prediction, or 



+ p\ + b- 



with a prediction of p= 1/2 or p=l . 
Where b and a are calculated by a least square fit to a linear curve given by 

iV /=-(JC-l)/2. 

and 

(K-D/2 

J] z * ^/+(/c-i)/2+r 

2 = izr (Ar_1)/2 



(K-l)/2 
i=~{K-\)l2 

for odd K. The indices are chosen so that the first Tj generated is multiplied by 
35 -(K-l)/2. 
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With reference to Figure 9, once the timing has been acquired, the next step is to 
apply a matched filter (step 96). The matched filter is applied using 1 sample per symbol at 
the filter output with no phase rotation. The next step is a second AGC acquisition on a block 
of 128 samples at 1 sample per symbol (step 98). 

At 1 sample per symbol, 32 samples are collected. The second AGC acquisition is 
performed on the block of 32 samples. During this step, the following is calculated. 

1 " 2 
Ems- — yizJ 

where N = 128 samples at 4 samples per symbol and the AGC value is calculated 
utilizing the following equation: 

Ems 



10 AGC = 



2xl6 2 



The next step is to perform fine frequency estimation (step 100). The fine frequency 
estimation is performed on the output of the matched filter after timing acquisition has been 
performed. Note that in order to increase performance at low SNRs, preferably two 
15 periodograms are averaged rather than using only one. A flow diagram illustrating the fine 
frequency estimation processing portion of the variable rate modem is shown in Figure 10. 

The following definitions apply: 

NFFT: the number of bins in the FFT, 256 for high SNR and 1024 for low SNR 
PADD: zero padding fraction, 0 for high SNR and 0.5 for low SNR 



With reference to Figure 10, the timing acquisition is performed utilizing the 
following pseudo code. 

Fine Frequency Acquisition Pseudo Code 



1 . for i= 0 to NFFT*PADD - 1 

collect the output of the matched filter (step 110) z / = /. : + jQ. at the correct 
time epoch (step 1 12); 

do z i - z* in two stages: 
30 calculate zf = n ; + jw n where w ( . = if - Qf, w t = 21 '.Q. 

calculate zf - x i + jy n where x t = uf - wf , y i = 2u i w i 

2. calculate a NFFT point of {z 4 with NFFT* PADD zeros padded (step 1 14). The 
absolute value of the FFT squared is denoted by ES. 

35 FFT shift ES: 
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ES'(i) = 



, NFFT~\ NFFT 
ES\ /- Li= ■» NFFT - 1 



W 2" 



NFFT 



J 



,i=0^ 1 

2 



where is the FFT shifted version of ES. From here on in, ES' is denoted by ES. 

( ) NFFT I 

3. find the maximum of {ES^.^ ' (step 116); 

4. find the index of the peak denoted by peak_index; 

5. the frequency estimation (step 1 18) is given as 

1 symbol rate ( f ■ . - , NFFT\ 
± = peak index 

4 NFFT \ 2 J 

6. alternatively, better results can be achieved if the modulation wipe off is performed 
as follows 



f 




4 j arctan 


(a)) 







exp„ t , _ Jf+Qf . ~ 

10 7. an approximation of the modulation wipe off can be achieved using the following: 



7.1 u=llog 2 (l? + Q?)\ 



This can be performed using iterated comparisons. The iteration is through the values 
-3 < v<4. 

15 7.2 The wipe off is performed by 

2 U 



With reference to Figure 9, the next step is to perform phase acquisition (step 102). 
The input to the phase acquisition portion of the method is Z n which are complex samples (I 
20 and Q) at a rate of 1 sample per second. The value of N is equal to 16 complex samples. One 

assumption made is that the frequency offsets are smaller than ^ ^ ^ where T is the 

symbol period, i.e., for N = 16, less than 1/128 of the symbol rate that is approximately the 
maximum frequency error. 

The phase acquisition process receives N complex samples (N = 16) and functions to 
25 determine the rotation angle that brings the samples after modulation wiping to a minimum 
variance. I and Q are the in phase and quadrature components, respectively. 

.z M =I H +jQ m9 n = 0.Jf-l u 
The method consists of two stages: (1) a coarse grained search that chooses between 4 
'hypothesis' of signal rotation with modulation unwiping in the first quadrant and (2) a fine 
30 estimation of the phase rotation angle. The role of the operation z 4 is to wipe off the effect of 
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the QPSK modulation on the phase. Before this operation, the phase of z is the sum of an 
unknown phase and of the phase due to modulation. After this operation, the phase depends 
less on the data. 

A flow diagram illustrating the coarse phase acquisition processing portion of the 
5 variable rate modem is shown in Figure 11. A flow diagram illustrating the fine phase 
acquisition processing portion of the variable rate modem is shown in Figure 12. 

With reference to Figures 11 and 12, the pseudo code for the phase acquisition 
follows. 

10 Phase Acquisition Pseudo Code 

1 . Coarse Grained Search 

1.1 for hypothesis Hk, k=l, 2, 3, 4,. ..rotate z n by adding the phase ©k = k-22.5° (step 
120) in accordance with 

15 z n (Q)= z.-e* 

1 .2 wipe off z n (00 modulation as follows (step 1 22): 

foreach z n =I n +jQ n 

check in which quadrant it appears 
if InQn < 0 then swap /„ and Q n 
20 take the absolute value of and Q n such that Z„ is now in the first 

quadrant; 

endfor 

1.3 sum the wiped off (or rotated) vectors z n (0$ using the following (step 124) (note 
25 that the sum should continue until N-l = 15) 

1.4 find the energy for each of the hypothesis (step 126) using 

E(® k )=\s(e k f 

At this point, it is relatively easy to show that for a rotational transformation, the 
30 minimum variance condition is equivalent to finding the maximum energy. 

1 .5 choose the hypotheses with maximum energy E(® kl ) (step 128). 

2. 'Fine' Phase Estimate 

2.1 let S(® k )= X + JY (step 130) be the vector having maximum energy; if Y > X 
3 5 (step 1 32) then swap X and Y (step 1 34) and set a swap Jlag (step 136); 

2.2 calculate a = atan(Y/X) using the techniques described below in 2.2.1 to 2.2.4 
(step 130); 

2.2.1 find a quadratic fit for the inverse term u=l/x: 
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'normalize' X to the range of [0.5. .1] which can be performed using left shift 
operations as follows 

x = 2 m -x 

a = 4.23493c 2 - 5.83 543c + 2.6192 

2.2.2 calculate 

v = ^ = 2--(y«) 

x 

where the 'unnormalization' by a factor of 2~ m is performed using right shifts (step 



137); 



2.2.3 calculate a = atan(v) using a quadratic fit (step 138): 
a = -0.0052v 2 + 1.0804v -0.2886 



10 Note: (1) these coefficients are scaled in radians; 

(2) this approximation assumes that 0 < v < 1 which is true because v=x/y and 
2. i assures that x > y; 

2.2.4 if the swapjlag is set (step 139) than a = 90° - a (step 140); 
2.3 the rotation angle 0, which must be subtracted from the phase of the samples z„, is 
15 given as 

® = ® k +a-45° (step 142) 



2.4 the rotation by © yields 



ze J 



This rotation operation brings in the received signal at phase ± it 1 4,±3^ / 4 , etc. in 
20 accordance with the data (for noiseless signals). Note that steps 132 through 140 comprise 
the detail implementation of step 130. 

The rotation angle 0 obtained during acquisition is the initial phase of the phase 
control loop. 

A block diagram illustrating pre-tracking and tracking portion of the variable rate 
25 modem is shown in Figure 13. This block diagram describes both the pre-tracking and 
tracking stages of the main loop processing. The I and Q data output from the matched filter 
is input to a mixer (multiplier) 211. The I and Q output of the mixer 21 1 is input to AGC 
estimation 213, timing detector 214 and phase detector 215. 

The AGC estimation 213 functions to generate an AGC control based on the energy 
30 or average power of the input signal. The AGC control is passed through a loop filter 216 
before being applied to the AGC hardware in the IF module section 1 6 (Figure 1). 

The timing detector 214 functions to generate a noisy timing estimate which is first 
filtered via loop filter 217 before being applied to the clock generation circuitry (not shown) 
in the A/D converter section of the baseband module 18. The A/D clock control signal output 
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from the loop filter 217 functions to adjust and correct for clock timing errors in the AID 
converters. 

The phase detector 215 functions to generate a noisy phase estimate that is input to 
the loop filter 218. One signal output of the loop filter comprises a VCO control signal which 
5 is input to the VCO circuitry (not shown) in the IF module 16. The VCO control signal 
functions to raise or lower the frequency output of the VCO. A second signal output of the 
loop filter 218 is the phase estimate <p that is fed into the e* (#) block 212. This block 
functions to correct for phase and small frequency offsets. The output of the block 
constitutes the second input to the I/Q mixer 211. 

10 The I and Q output of the mixer 21 1 is also input to the soft decision block 219. The 

matched filter outputs z n , after time synchronization, frequency and phase correction, must 
pass the soft decision block 219 before being input to Viterbi detection in the Viterbi decoder 
(not shown) in the baseband module 18 (Figure 1) (see also step 158 Figure 14). This process 
is performed on I n and Q n independently. The Viterbi decoder receives for each I n and Q n , 3 

15 or 4 bits of quantization. The exact number depends on the specific Viterbi detector. One of 
the bits is the sign and the other bits describe the signal level. 

For example, if the level thresholds are cti < a.2 < (X3, the two bits of soft decision 
which describe the level are shown below in Table 7: 



Table 1 


In 


Bit Representation 


In < CCl 


00 


ai < I„ < a 2 


01 


a 2 < I n < a 3 


10 


a 3 < In 


1 1 



With reference to Figure 3, once the acquisition stage is complete, pre-tracking can be 
performed (step 37). The pre-track stage is in principle normal tracking of AGC and timing, 
while the phase tracking loop has a larger bandwidth and the 90° ambiguity still must be 
25 resolved. A flow diagram illustrating the pre-track processing portion of the variable rate 
modem is shown in Figure 14. 

The method comprises (1) AGC tracking as described hereinabove, (2) timing 
tracking using a method that operates on a group of symbols (at one sample per symbol) to 
generate an estimate that is filtered using a single pole IIR filter; the method is called at a rate 
30 of approximately 160 Hz which 32 symbols being processed each time, and (3) phase 
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tracking performed once per symbol for symbol rates below 64 KHz and once every 2 or 4 
symbols for higher symbol rates. The criteria for declaring tracking lock can be a timeout, a 
high raw BER or one or more lock detectors. 

First, the l and Q data is input to the matched filter (step 150). AGC tracking is then 
5 performed on the output of the matched filter (step 152). A flow diagram illustrating the 
AGC tracking portion of the variable rate modem is shown in Figure 15. The pseudo code 
for the AGC tracking follows. 

AGC Tracking Pseudo Code 

10 1. initialize variables including 'count' and set the threshold to 256 (step 180); 

2. check if count = 8 (step 182); 

3. every T milliseconds (T « 100 ms) receive the output z n {n represents one sample 
per symbol) of the matched filter (step 186) while count is < 8; 

4. calculate g„ as follows (step 196) 

15 g n = Treshold - \z n \ 2 

5. calculate AGC value (step 198). Filter g n using a Butterworth low pass filter (LPF) 
transfer function shown below 

0.073 •(! 4- z' 1 ) 



H(z) = 



1- 0.854- z* 1 



with a 3dB cutoff rate at 0.1 of the digital Butterworth LPF. The corresponding 
20 difference equation is as follows 

G nH = 0.854 • G„ + 0.073 • (g n + g^) 

6. the output Y„ of this filter is the new 'inverse AGC (in the energy domain) (step 

200). 

AGC = — 

25 7. to prevent clipping at the input to the A/D converter, the threshold is adapted every 

K*Tms(K«8, 16). 

8. find Y n in an inverse _AGC_table and write the value to the AGC circuitry in the IF 
module; 

8. increment count by one; - ■ " : 
30 9. after N samples (N = 8 in the above example) reset count to 0 and calculate Ems 

(step 184); 

10. if Ems > 0.25 (step 188) then increase the threshold (step 192); 

Threshold = a Threshold, a ^ 1 + 1/8 
35 else if the threshold > 1/256 AND Ems < 0.24 (step 1 90) then decrease threshold 

(step 194); 

Threshold = /J - Threshold, (5 ^1-1/32 . 

endif 
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With reference to Figure 14, the next step is to perform timing tracking (step 154). 
The timing tracking portion comprises a loop filter, an ideal integrator and an inherent 
numerically controlled oscillator (NCO): A block diagram illustrating a physical model of 
5 the timing tracking portion of the variable rate modem is shown in Figure 16. The filtered I 
and Q data at 1 sample per symbol, i.e., the output of the mixer 211 (Figure 13), is input to a 
Mueller error detector 170. The output of the error detector 170 is input to a loop filter 172. 
The output of the filter constitutes the A/D clock control signal. This signal is input to NCO 
174 that functions to generate the clock for the A/D converter. The output of the NCO is also 
10 fed back to the Mueller error detector 1 70. 

The timing signal model of the loop filter 172 can be represented as the following 



l-z 



and the NCO can be represented as 



1 

l-z' 1 



15 Note that the loop filter has a gain of p, an ideal integrator and a zero at y, i.e., proportional 
and derivative control. The timing tracking loop utilizes the Mueller and Mueller timing 
error detector 170 which operates on the filtered in-phase and quadrature components on a 
block of 32 symbols at a rate of approximately 320 Hz (determined by the operating system). 
The timing error detector is given by the average of the I and Q error components: 



20 e = 



2K 
where 

€m = In • sgn(/ fl _, )-/.-,• sgn(/ n ) 
£ Qn = Qn • sgn^)- fi,-t • sgn(0„ ) 
and K - 1/16 is the nominal value of the signal amplitude as set by the AGC. Next, the 
25 average of a block size N = 32 is computed as follows: 

1 Af-l 

The error process is filtered by the loop filter 172 as follows: 

t n =t n _, + f}{E n -yE n _ x ) 
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where t n is the correction written to the hardware. The NCO 174 or ideal integrator is part of 
the system and may or may not be implemented in software. It should, however, be taken 
into account when calculating the closed loop transfer function. Based on the timing signal 
model given above and the assumption that the error detector is ideal, i.e., it yields the exact 
timing error), the closed loop transfer function of the estimated timing epoch is given by 

T l + (/3-2)z- 1 +(l-/3y)z- 2 
Note that it is convenient to construct the filter using the classical parameters from 
continuous time in terms of co n = 27if n , the natural frequency and the damping factor, it is 

desirable to express p and y in terms of co n and This can be achieved in two stages. The 
first stage is to convert to the discrete time domain standard form by solving the following set 
of equations: 

\{J3 - 2)z~ x =" -2r cos Qz~ l 
i(l-^>- 2 =rV 2 



from which we obtain 



9 = cos" 



2r 



15 or for the other direction 



P = 2(\-rcos6) 

1-r 2 1-r 2 

7 = 



P 2(l-rcos<9) 

We now express r and #in terms of the design parameters a> n and § as follows. 



20 This is described in more detail in Best, Phased Locked Loops, Second Edition, p. 342, 
McGraw Hill, 1993. In addition, a more detailed description of the timing error detector 170 
is given in IEEE Transactions on Communications, May 1976, pp. 516 to 530. 

The following is based on the assumption that the timing error detector gain is unity. 
The detector gain will now be discussed in more detail. The transformation is based on the 

25 impulse invariance conversion from continuous to discrete time. From r and 6 we can 
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compute y?and y Based on the above equations, assuming an operating rate of 1/T = 320 or 
640 operations per second and setting £ = 0.707, the parameters for J3 and y are presented 
below. 

5 Table 8: Parameters /? and y 



f n (Hz) @ 320 Hz 


P 


Y 


/„ (Hz) @ 640 Hz 


0.1 


0.0028 


0.9986 


0.2 


0.2 


0.0056 


0.9972 


0.4 


0.32 


0.0089 


0.9956 


0.64 


0.5 


0.0139 


0.9931 


1 


1 


0.0278 


0.9863 


2 


1.6 


0.0444 


0.9783 


3.2 


2 


0.0555 


0.9730 


4 


3.2 


0.0888 


0.9575 


6.4 


5 


0.1386 


0.9351 


10 


7.5 


0.2075 


0.9058 


15 


10 


0.2760 


0.8784 


20 



The choice of £ = 0.707 is due to the dependence of the gain on the error detector which is a 
function of the AGC and the signal to noise ratio. 

With reference to Figure 14, the next step after timing tracking is phase tracking (step 
10 156). A block diagram illustrating the phase tracking portion of the variable rate modem is 
shown in Figure 17. The phase tracking comprises three main components including a phase 
detector 254, loop filter 258 and a frequency jump 262 if the frequency is larger than a 
threshold. 

The phase tracker mainly operates at the symbol rate. For higher data rates, real time 
15 constraints may pose a problem. For high symbol rates it may be necessary to operate the 
phase tracking at rates slower than the symbol rate. 

The filtered I and Q data is input to a mixer 252. The output of the mixer 254 is input 
to the phase detector 254. The QPSK phase detector performs the following: 

e n =Q n sgn(I n )-I n *sgn(Q n ) 
20 where I and Q denote in-phase and quadrature components, respectively, at a rate of once per . 
symbol. 

The output of the phase detector is input to a loop filter 258. The filter is presented in 
state space representation. The filter is defined by the parameters ao and a;. The Doppler 
estimate, D n , is preferably also input to the filter from an external loop. The equations are 
25 shown below wherein F n and P„ denote the frequency arid phase, respectively. Note that the 
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phase P„ is the same phase discussed earlier denoted <p or 0. The filter preferably comprises 
an infinite impulse response (IIR) filter. The loop filter performs the following: 

F^^a l {P n -P n _ l ) + b l F n ^D n 
P n «=b Q e n +P n +F R+i -D n 

where s n is the phase error from the phase detector and ai = I - bi. 
5 If the frequency F n in the filter is large, a jump in the direct digital synthesizer (DDS) 

or other similar frequency source is performed having a size Af. A typical value of Af 
corresponds to 10 Hz. 

The absolute value of the output of the loop filter F n+I is compared against a threshold 
via a compare block 260 that performs the following comparison: 

io K.l>f 

If the expression is true then jump Af in the opposite direction via frequency jump 262. Note 
that the F n and threshold values are in units of Radians/operating period. 

The data is written to the direct digital synthesizer (DDS) using frequency jump unit 
262 in accordance with the following 
1 5 synthesizer _ freq — synthesizer _ freq + Af 

Note that the quantity Afis in units of Radians/operating period. 
The following variables apply: 

P^P.-^sgnm 
F n =F n -Af-sgn(F n ) 

It is assumed that the Doppler variable is zero. For a non-zero Doppler value one skilled in 
20 the art can modify the above equations to include the Doppler variable. 

Note also that the Doppler estimates are supplied from an external loop uncoupled to 
the IIR filter. 

The output of the frequency jump 262 and the compare 260 are input to the e? (m) block 
264. The output of block 264 fonns the second input of the mixer 252. 
25 The transfer function of the IIR used in the loop filter 258 will now be described in 

more detail, by taking the Z transform of the two IIR equations (F' n +h Pn+i) and substituting 

e = '6 -P 

where <j> n is the phase of the input signal (linear approximation), we obtain: 
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0. _ b 0 {\-b x z- x ) 

P n l-(l + a 0 )z- 1 +(a 0 6 l +« I )z- 2 

Note that it is convenient to use the normalized form for a second order system, i.e., to define 

V2 

the system by £ = 1 or and co n which is determined by the required lock range, Doppler 

offset or output signal to noise ratio (SNR). A damping coefficient of £ = 1 is recommended 
. 5 for reducing cycle slips as described in Ascheid and Meyr, "Cycle Clips in PLLs: A Tutorial 
Survey," IEEE Transactions in Communications, Vol. 30, October 1982, pp. 2228 to 2241. 
For acquisition it is preferable to use a smaller damping coefficient in the range 

V2 

£ = 0.5 to which yields a faster response with an overshoot of 1 8% for £ = 0.5 . 

After phase tracking is accomplished, the I and Q soft decisions are output to the 
.10 Viterbi decoder (not shown) (step 158). Whether signal lock has been achieved is then 
determined by comparing the BER to a predetermined threshold (step 160). If lock has been 
achieved, the main tracking loop is entered. If it is determined that lock has not been 
achieved than the system goes back to the acquisition (step 161). 

A flow diagram illustrating the tracking processing portion of the variable rate modem 
15 is shown in Figure 18. The various functions performed in tracking are substantially the 
same as those performed during the pre-tracking stage except for the phase tracking which 
has a smaller loop bandwidth. 

The tracking stage comprises the steps of performing a matched filter (step 220), 
AGC tracking (step 222), timing tracking (step 224), phase tracking (step 226), I and Q soft 
20 decisions (step 228) and determining whether lock has been achieved (step 230). If it is 
determined that lock is achieved, the signal detection is declared (step 232). 

In operation, the main loop advances from stage to stage until tracking is achieved. 
At that point, the BER is measured by the Viterbi decoder. If the BER indicates failure, then 
the system returns back to acquisition. 
25 If the BER indicates lock, then the system remains in tracking until the BER indicates 

that the system is not in lock due to a large number of errors. Other indications may be used 
such as the behavior of phase locking loops. 

While the invention has been described with respect to a limited number of 
embodiments, it will be appreciated that many variations, modifications and other 
30 applications of the invention may be made. 
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